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Abstract 

\ In this talk it is discussed the derivation of low-frequencies part of quark determinant 
and partition function. As a first application, quark condensate is calculated beyond 
chiral limit with the account of 0(m), O(^), O(jj-m) and 0(jj-m In m) corrections. 
It was demonstrated complete correspondence of the results to chiral perturbation 
theory. 
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2 : Introduction 

Instanton vacuum model assume that QCD vacuum is filled not only by perturbative 
but also very strong non-perturbative fluctuations - instantons. This model provides a 
natural mechanism for the spontaneous breaking of chiral symmetry (SBCS) due to the 
derealization of single-instanton quark zero modes in the instanton medium. The model 
is described by two main parameters - the average instanton size p ~ 0.3 fm and average 
inter-instanton distance R ~ 1 fm. These values was found phenomenologically |TJ and 
theoretically [2] and was confirmed by lattice measurements jSJ EJ |5J EH U\ ■ On the base of 
this model was developed effective action approach 0121110], providing reliable method 
of the calculations of the observables in hadron physics at least in chiral limit. 

On the other hand, chiral perturbation theory makes a theoretical framework incorpo- 
rating the constraints on low-energy behavior of various observables based on the general 
principles of chiral symmetry and quantum field theory [llj . 

It is natural expect, that instanton vacuum model leads to the results compatible with 
chiral perturbation theory. 

One of the most important quantities related with SBCS is the vacuum quark con- 
densate < qq >, playing also important phenomenological role in various applications of 
QCD sum rule approach. Previous investigations 12J shows that beyond chiral limit and 
at small current quark mass m ~ few MeV these quantity receive large so called chiral 
log contribution ~ J-mlnm with fixed model independent coefficient. On the typical 
scale lGeV it become leading correction since l-r^ln ml > 1. It was shown, that this 

^^c/y c 

correction is due to pion loop contribution [T2] ITT] . 

So, to be consistent we have to calculate simultaneously all of the corrections of order 
m, jj^, jjr- In m in order to find quark condensate beyond chiral limit. 



1 Talk given at the XVII International Baldin Seminar " Relativistic Nuclear Physics and Quantum 
Chrome-dynamics" , Sept.27-Oct.2, 2004 (JINR, Dubna, Russia). 



1 



In our previous papers O Ej on the base of low-frequencies part of light quark 
determinant Deti ow , obtained in ^5J|B1EI|, was derived effective action. In this framework 
was investigated current quark mass m dependence of the quark condensate, but without 
meson loop contribution [T^j . 

In the present work we refine the derivation of the low-frequencies part of light quark 
determinant ~Deti ow . The following averaging of Deti ow over instanton collective coor- 
dinates is done independently over each instanton thanks to small packing parameter 
7r(^) 4 ~ 0.1 and also by introducing constituent quarks degree of freedoms ip. This proce- 
dure leads to the light quarks partition function Z[m]. We apply bosonisation procedure 
to Z[m], which is exact one for our case Nf = 2 and calculate partition function Z[m] 
with account of meson loops. This one provide us the quark condensate with desired 
0{m), 0{jj-), 0(jj-m In m) corrections. 

Low— frequencies part of light quark determinant 

The main assumption of previous works [HI EH EH (see also review ^H]) was that at very 
small m the quark propagator in the single instanton field Aj can be approximated as: 

c t n\ 1 , l $ 0/ >< $0/1 m 
SAm ~ 0) « — + 1 (1) 

id »m 

It gives proper value for the < &oi\Si(m ~ 0)|$o/ >= but in S/(m ~ 0)|$oj >= 
+ ^l$o/ > second extra term has a wrong chiral properties. We may neglect by this 
one only for the m ~ 0. 

At the present case of non- small m we assume: 

Si « S + 3 id >< ^ idS , S = (2) 

cj id + im 

where 

Cl — — < $ 0/ |ii95oi<9|<I>o_f >= im < <& OI \S id\$oi > (3) 
The matrix element < &oi\Si\&oi >— — , more over 

1 1 

5/ $ M >= — $o/ >, < $oi\Si =< $o/|— (4) 
im im 

as it must be. 

In the field of instanton ensemble, represented by A = J2i Ai, full quark propagator, 
expanded with respect to a single instanton, and with account Eq. (J2J) is: 

S = So + - So) + Y,( S i- So)S \Sj - S ) 

I I^J 

+ E (Si ~ So)S x (Sj — S )S x (Sk — So) + ... 
= So + Y, Soid\%! > (i + ^ri + ...)_< $ j\idS 



i, j 



i j 



S + Y,S id\% I > (c^r) /J <$oj| ^ (5) 



IJ 
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where 

Cu = Sjjci = -S u < <$> 0I \idS id\ ® 0I >, 

(C - T)jj = -< $ 0I \idS id\$ j > (6) 

We are calculating Deti ow using the formula: 

lnDet/ ^ = Tr / idm'(S(m!) — So(m')) (7) 

J Mi 

Within zero-mode assumption (Eq. (J2J)) the trace is restricted to the subspace of instan- 
tons: 

Tr(S - So) = -£ < <S> 0; j\idS 2 id\<f> 0l i >< $o,/|(-^)|<£o,j > (8) 
i, j idS id 

Introducing now the matrix 

B(m)u =< * 0tI \idS o id\$ ,j > (9) 

it is easy to show that 



rm rB(n~i) j 

lnDetz ^ = Tr / idm'(S(m) - S (m')) = V / (dB(m) 

JMi jJb(Mi) B{m!) 



in 



B (tti) 

= Trln = lndetS(m) -lndetS(Mi) (10) 

which is desired answer. The determinant det B(m) from Eq. (jlUj) is the extension of the 
Lee-Bardeen result |T5] for the non-small values of current quark mass m. 



Light quark effective action beyond chiral limit 

Averaged Deti ow leads to the partition function Z[m], which for Nf = 2 has the form: 

2 

Z[m] = f dX+dX-Di/jD^ exp[f d*xJ2^ f (id + im f )ip f (11) 

+X+Y+ + \_Y~ + N + In — + AL In — ], 

A_ 

here X± are dynamical couplings (K is unessential constant, which provide under- logarithm 
expression dimensionless) jHJIISl El- Values of them are defined by saddle-point calcula- 
tions. Y% are t'Hooft type interaction terms |lUj : 

Yf = — — / d 4 x\(l - —) det iJ^ip, x) + —— det iJ±(x)] (12) 
2 N^-lJ LV 2N C J J 8N C ^ n v ; 
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where q(k) = 2irpF(k)ip(k). The form-factor F(k) is due to zero-modes and has explicit 



x 



form F(k) = - ±[I (t)K (t) - h(t)Kx(t)] t=M p. In the following we will neglect by J± v fg 
interaction term, since it give a O(^) contribution to the quark condensate. Since 
Q( x ) = I exp(ikx) q(k), jf g (x) = qj(x)^q g (x), 

, i.J + (x) . iJ~(x) ,„„,. 

det — + det — (13) 

9 9 

= ^(~(<l + ( x )<l( x )) 2 - (q + (x)i<y 5 Tq(x)) 2 + (g 4 (x)rq(x)) 2 + (q + ' (x)i-f 5 q(x)) 2 ) . 
8g 2 

Here color factor g 2 = ^^iv-i^ • 

In the following we will take equal number of instantons and antiinstantons iV + = 
iV_ = N/2 and corresponding couplings A± = A. 

Now it is natural to bosonize quark-quark interaction terms (113(1 by introducing meson 
fields. For Nf = 2 case it is exact procedure. We have to take into account the changes 
of q and q^ under the £77(2) chiral transformations: 

Sq = ijsraq, 5q + = q + i^Ta 

to introduce appropriate meson fields, changing under 577(2) chiral transformations as: 

5a = 2a<f), 5(f) = —2aa, 5n = —2aa, 5a = 2i]d. 

Then 5q + (a + ij 5 f(j))q = 0, 5q + {fa + ijs^q = means that these combinations of fields 
are chiral invariant 2 . So, the interaction term has an exact bosonized representation: 



/ d 4 xexp[A(det — + det — )] (14) 
J a a 

A^ 
2g 



r - t A ' 5 - 1 

J DaD(f)DT]Daexp J d 4 x[-^—q + i(a + %^^f(ft + ira + ^rr)q - -(a 2 + (ft 2 + a 2 + if)} 



Then the partition function is 

K 



Z[m] = J d\DaD0Dr]Da exp[N In — - N (15) 
J d 4 x(a 2 + ft 2 + a 2 + r] 2 ) + Tr In 



1 f * ( 2 ^ 72 ^_ ^2 j_ 2\ 1 rj-L, 1 P + %m + ^( 27r P) 2F ( ff + ^0 + ITa + l5 T])F 

- / a x(a + d> + a +77 ) + Tr In 

2 J y 1 ' p + im 1 



(Tr(...) means here tr 7jC j / d 4 x < x\(...)\x >, where tr 7iC j is the trace over Dirac, color, 
and flavor indexes.) In the following we assume m u = = m. Then common saddle 
point on A, a (= const) (others = 0)is defined by Eqs. = dv[mx<r} _ w Y iere 

the potential 

V[m,X,a] = -AT h, g + N + V - IV In t± ''<'" ±M£l <W> (16) 
A 2 p + im 



2 Certainly, quark-quark interaction term Eq. (|13fl is non-invariant over U(l) axial transformations, 
as it must be. 
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and we denned M(X, a) = ^-{2-Kpfa. Then the common saddle-point on A and a is 
given by Eqs.: 

JV = iTV = i Vo .. (17) 

2 p + i(m + M(A,cr)F 2 (p)) 2 V ; 

The solutions of this Eqs. are A and cr = (2-j^) 1 / 2 = 2 1 / 2 i?~ 2 . It is clear that M = 
M(Ao,cr ) has a meaning of dynamical quark mass, which is defined by this Eqs.. At 
typical values Br 1 = 200 MeV, p~ l = 600 MeV we have a\ = 2(200 MeV) 4 , and in chiral 
limit m = M — > M 00 = 358 MeV, Aoo ~ Mq . It is clear that due to saddle-point 
equation (fTTj) M (and Ao) become the function of the current mass m. This dependence 
was investigated in [14J. 



Vacuum with account of quantum corrections 

The account of the quantum fluctuations around saddle-points cr ,A will change the 
potential V[m, X, a] to V e ff[m, X, a] (it is clear that the difference between these two 
potentials is order of 1/N C ). Then, the partition function is given by Eq. 



Z[m] — J dXexp(—Veff[m,X,a]) (18) 

There is important difference between this instanton generated partition function Z[m] 
and traditional NJL-type models - we have to integrate over the coupling A here. As was 
mentioned before, this integration on A by saddle-point method leads to exact answer. 
This saddle-point is defined by Eq.: 

dVpff \m, X, a] , , 

etn ' 1 = (19) 
dX 

which leads to the A as a function of a, i.e. A = A(cr). 

Then, the vacuum is the minimum of the effective potential V e ff[m, a], which is given 
by a solution of the equation 

dVeff[m,a,X(a)] = dV eff [m, a, A (a)] = Q 
da da 

where it was used Eq. (IT9*|) . 

We denote a fluctuations as a primed fields The action and corresponding V e ff 
now has a form: 

S[m, A, a, $'] = S [m, A, a] + S v [m, A, a, $'], (21) 

or ^ 1 T rr ^ 1 1 t t 2 rrn i P + Um + M (X, a) F 2 ) K Ar 

S [m,X,a] = V\m,X,a] = -Va 2 - Trln- - — — - - iVln — + N 

2 p + im A 



S v [m, X, a, $'] = J d 4 x^(a' 2 + $' 2 + a' 2 + V ' 2 ) 



1 iM(X,a)F 2 -a,—*, 

rlr r^r{a + lJ 5 T<p + IT a + 



2a< 



and 



p + i(m + M(X,a)F 2 ) 

V eff [m, X, a] = S [m, X, a] + V^ s \m, X, a] (22) 
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Here second term in Eq. (|22j) is explicitly represented by 

d 4 p 



V™?\m X a] - i T r In ^y[m,X,a^] _V r tfq_ _ J_ r 



M(A,a)F 2 (p) M(A,tx)F 2 (p + g) 

p + i(m + M(A,(x)F 2 (j9)) J j3 + g + i(m + M(A,(x)F 2 (p + g)) * J ' 1 J 

where the factors Tj = (1, 275T, ir, 75) and the sum on i is counted all corresponding meson 
fluctuations a',<fi',a',r}'. tr here means the trace over flavor, color and Dirac indexes. 
Integrals in Eq. (|23|) are completely convergent one due to the presence of the form- 
factors F. 

Certainly the quantum fluctuations contribution will move the the coupling A from Ao 
to A + Ai and a as cr ~* °"o + where ^ and ^ are of order 1/N C . 
First, consider Eq. ()19|): 

dV e// [m, A, gj _ 1 zM(A, a)F 2 V f A 

cU 2 p + i(m + M(A,a)F 2 ) 2 (2?r) 4 1 J 

2 + f d A p M(X,a)F 2 (p) r M(X,a)F 2 (p + q) 



x a — tr / v y — r r 

1 i (27r) 4 p + i(m + M(A,a)F 2 (p)) *p + g + i(m + M(A, a)F 2 (p + q)) 

d 4 p M(X,a)F 2 (p) ^ M(X,a)F 2 (p + q) 



x[-tr J 



(2tt) 4 p + i(m + M(A, a)F 2 (p)) l p + q + i(m + M(A, (j)F 2 (p + g))" 



+zt r f ^ ( M(X,a)F 2 (p) \ 2 p M(X ) a)F 2 (p + q) 

1 1 7 (2vr) 4 \p + i(m + M(X,(T)F 2 (p))J 1 p + q + i(m + M(X,a)F 2 {p + q)) tl 

From this saddle-point Eq. we get A = A(cr). 

From vacuum Eq. (J2U|) we in similar manner arrive to: 

dV eff [m,a,X(a)] 2 iM(A(cr), o~)F 2 V /• rf 4 g 

da ' p + z(m + M(A(a),a)F 2 ) (2i) 4 1 ' 

xf 2 _fr / ^ P M(A(a),a)F 2 (p) M(A(cr), a)F 2 (p + g) 

L 7 (27r) 4 p + i(m + M(A(a),a)F 2 (p)) l p + q + i(m + M(A(cr), a)F 2 (p + q)) A 

x[2 ? tr ( M{X{ala)F 2 {p) \\ M(X(a), a)F 2 (p + q) 

J (2n) 4: \p + i(m + M(X(a), a)F 2 (p)) J ' p + q + i(m + M(X(a), a)F 2 (p + q)) A 

Since we are believing to jj- expansion, it is natural inside quantum fluctuations contri- 
bution (under the integrals over q) to take a = a , M(X(a), a) = Mq. 

To simplify the expressions introduce vertices V^g), V$i(q) and meson propagators 
Ilj(g), which are defined as: 

V2M) = tr / ; , Ti 26 

J (2tt) 4 p + i/io(p) p + g + «/i (p + g) 

V + r, (27) 

J (2vr) 4 \p + ifJ>o{P)J p + q + %Ho{p + q) 

nr\q) = ^-V 2l (q). (28) 
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Here M (p) = M F 2 (p), Ho(p) — tti + Mq(p) and was taken into account that a 2 = 2R 4 . 
From Eqs. (|24j) and (J25j) we have 



Mi 
Mo 



2 1 / M (p) 
— — r + — Tr ' 



i? 4 1/ Vp + ^oCp), 



E / t^WC^-Vk^^C?) (29) 



- = -?E/7^^(g)n,(g) (30) 



a 4 (2tt 

The vertices V2i(g), V^g) and the meson propagators Hi(q) are well defined functions, 
providing well convergence of the integrals in Eqs. (jZSj) ) (jHO))- 

It is of special attention to the contribution of pion fluctuations ft at small 
pion momentum q. We shall demonstrate that this contribution leads to the famous 
chiral log term with model independent coefficient in the correspondence with previous 
calculations in NJL- model [TS] . 

Pion inverse propagator of U.^(q) at small q ~ m n is: U.^ t (q) = f^ in (m 2 + q 2 ). At 

lowest order on m, fkin,m=o ~ fn = 93 MeV, m 2 ~ m. 

The vertices in the right side of Eq. (j2§j) at g = and in chiral limit are: 

2 r d 4 v v 2 M 2 (v) 

V^{0) = ^I, rti^O) - ^;,^„(0) = 8^/ (^I (/ + m^ ))2 (31) 



We see that the factor in the left side of Eq. (|29j) in the chiral limit is equal to: 

d^p ipM (p) 



tr 



/ <^CT^ - "^ W» " W>» (32) 



(2n)*(p + tM (p)y 

Collecting all the factors we get small q < k contribution of pion fluctuations ft: 
a lt Mi, 3 /■« d 4 q 1 



aQ \4>>, small q ' Mq \#, small q ' 2 J2 J Q (2tt) 4 77^ + g 2 



(33) 



3 f K 2 j 2 ' 3 , r, , 771 2 



1 ^ -JT7 o\ = ~ oo 2 12 ^ +^7r m ' 



327T 2 f 2 Jo * * f 2 (ml + q 2 ) 327r 2 / 2 v 77 k 2 + m 2 



Here we put m = everywhere except m f . We see that the coefficient in the front of of 
m 2 In m 2 is a model independent as it must be. 
Quick estimate, assuming k = p _1 , gives 

ai < - Ml < -■ 3 ■ _2„2i„™2„2 



C7 ^' Mo 1 ^' ' 32tt 2 /> 2 



1+myinmy) « -0.4(1 + 0.054 In 0.054) (34) 



So, we expect that pion loops is provided not only non-analytical -^m In m term but also 
very large contribution to jr term. 

This one dictate the strategy of the following calculations of <j\ and M\. 

1. we have to extract analytically jmlnm term from pion loops; 

2. rest part of o\ and M\ can be calculated numerically and expanded over m, paying 
special attention to the pion loops and keeping and jfWn terms. 
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For actual numerical calculations we are using simplified version of the form-factor 
F(p) from ^7] (with corrected high momentum dependence): 



F(p < 2GeV) = — ^-j, F(p > 2GeV) = ^ (35) 
L z + p z p 6 



where L^^ = 8A8MeV. 



p 

At N c = 3 semi-numerical calculations of Mi and o\ lead to: 

M 

— = -0.662 -4.64m -4.01m lnm (36) 
M V ; 

— = -0.523 -4.26 m -4.00 m lnm (37) 

Here m is given in GeV . Certainly, in ()36|) the m lnm term is completely correspond to 
Eq. ()33|). is —66% in chiral limit and reach its maximum ~ —20% at m ~ 0.115 GeV. 

The relative shift of the vacuum — is —52% at the chiral limit and reach its maximum 
~ -2% at m ~ 0.125 GeV. 

The main contribution to both quantities ^ and ^ come from pion loops. Other 
mesons give the contribution ~ 10% to 0(^-) and O(jj-m) terms. 



Quark condensate 

We have to calculate quark condensate beyond chiral limit taking into account 0(m), 
O(-jkr-), O(-^-m) and O(-^mlnm) terms. Quark condensate is extracted from the partition 
function: 

1 dV ef f[m,X,a] 1 d(V[m,X,a]+V^f[m,Xo,a ]) 

< an > = — = — 

2V dm 2V dm 

= -±Tr(^-L- - ^!_) + _L av ?/'[m, Ao, „] 
2V p + ip,yp) p + im 2V dm 

here A = A + Ai, , a = a + a u M = M + M u /x(p) = m + MF 2 (p). First term of 
Eq. (jUD is 

~Tr{ - 1 -. r - -^) (39) 
2K p + i/i[p) p + im 

AN f ^ ( ^ m M 1 M 0(P)(P 2 -^0(P)) \ 

c 7 (2vr) 4 l^p 2 + ^g(p) p 2 + m 2 + M (p 2 + /i 2 (p)) 2 J 

Second term of Eq. (J3%|) - meson loops contribution to the condensate is 

1 gK?f[m,A ,a ] i r d*q ( [ d*p M (p) M (p + g) \ 

2V 9m 2 YM2vr) 4 V V (27r) 4 (p + i/i (p)) 2 l £ + g + z>o(p + q) J 

y (2vr) 4 p + ^ob) P + g + i>o(p + g) 7 v ; 



s 



At m = and without meson loops the condensate is 

f d A p M 00 (p) , s 

Here M 00 = M 0)m=0 . 

Let us to consider now the contribution of pion fluctuations </>' to the quark condensate 
at small q. First we consider: 

1 dV e f/™ H «[m, A , a ] _ r d*p M*(p)M f d\ 



9m c 7 (2tt)Mp 2 + ^(p)) 2 Vo (2tt)VL« + <7 2 ) 

We keep m only in m. 2 . Then at m = fMoip) =>■ =>■ Mx)(p), /fcin =>■ /tt and we have 

<qq> = <qq> 00 -—AN c J (2?r)4 (p2 + Mo2()(p))2 (43) 



+ 12 AT, 



d 4 p M 3 (p) /•« d 4 g 1 



(2vr) 4 (p 2 + M 2 (p)) 2 7 (2vr)4 / 2 (m 2 + q 



-,2^ 



A 3 r A i \ 

=< qq > 00 1 - - / , u , 2 , 2 - 2 x 44 
V 2 Jo (2tt) 4 f£(mi + q 2 )J 

Eq. ()33j) for -jg 1 was applied here. We see that Eq. (|4^|) is in the full correspondence with 

mum 

Detailed numerical calculations lead to the semi-analytical formula for the quark con- 
densate including all 0(m), O(jj-) and 0(j=j-m In m)-corrections: 

< qq >=< qq > m= o (1 — 18.53 m — 7. 72 m In m) (45) 

Here < qq > m =o= 0.52 < qq > 00 . Certainly, the mlnm term in Eq.(|43j) is in full 

correspondence with Eq. (JUJ), as it must be. < qq > / < qq > m =o is a rising function of 

m until m ~ 0.04 GeV and is a falling one in the region m > 0.04 GeV. 

The main contribution to 0(4r), O(-^m) and Of^-mlnm) terms in J^p" is due to 

\ N c n ^N c > \ N c ' <qq>oo 

pion loops. Other mesons give the contribution ~ few% to O(jj-) and O(-^-m) terms. 

md — rn u effects in quark condensate 

Current quark mass become diagonal 2x2 matrix with mi = m u ,m 2 = m^, m = 
mi^y 1 + = m + Sm 1 ^. Here m = mi + m2 ; = mi — m 2 . Let us introduce 

external field s«. In our particular case it is S3 = z^p, si = s 2 = 0. Our aim is to find the 
asymmetry of the quark condensate <uu ^^ d> , taking into account only 0(5m) terms 
and neglecting by 0(^-5m), O(^-Smlnm). It means that we neglect at all by meson 
loops contribution. 

In the presence of the external field s we expect also vacuum field o. Effective potential 
within requested accuracy is 

V eff [a, a,m) « S [m, A, a, a] (46) 

= V + ^ _ Tlh ptift-ti(-tf(V:.y) _ + iv. 

2 p + im + zrs A 



\,a,a are defined by the vacuum equations: 



9V eff n dV eff dV eff 

^r = ' ^ = ' ^r = - (47) 

They can be reduced to the following form: 

l Ti F 2 (p)M t (m l + M l F\p)) = 

2 p 2 + (m, L + MiF 2 (p)) 2 { ' 

where Mi = ^-(2np) 2 (a ± 0-3). Solution of these equations leads to A = X[m, s\, a = 
a[m,s\ Oi = ai[m,s\. We have to put them into V e ff and find V e ff = V e ff[m,s\. Desired 
correlator is 

dV e ff[m,s\ 

ds3 L 3 =^, Sl , 2 =o (49) 

We calculate this correlator within requested accuracy, taking into account only 0(5m) 
terms. So, the difference of the vacuum quark condensates of u and d quarks is 



<uu> - < dd> (50) 
1 

~ V 



Tr(^ — — — - - : ) - Tr( 



K p + i(m u + M U F 2 ) p + im u p + i(rrid + M^F 2 ) p + irrid' 

We expect that < dd > << uu > if > m u . 

Typical values of light current quark masses JH] are m u = 5.1MeV, = 9.3MeV on 
the scale lGeV (which is in fact close to our scale p _1 = 0.6 GeV) leads to the asymmetry 

= 0.026 (51) 

< uu > 

From this asymmetry and using sum-rules [20] we estimate strange quark condensate at 
m s = 120 MeV as: 

= 0.43, (52) 

< uu > 

which is rather small. The reason that the asymmetry ()51|) is rather large. 



Conclusion 

In the framework of instanton vacuum model it was calculated simplest possible correlator 
- quark condensate with complete account of 0(m), 0(4-), 0(4-m) and O(-^-mlnm) 
terms, demanding the calculation of meson loops contribution. Since initial instanton 
generated quark-quark interactions are nonlocal and contain corresponding form-factor 
induced by quark zero-mode, these loops correspond completely convergent integrals. 
The main loop corrections come from the pions, as it was expected. We found that 
O(jj-) corrections are very large ~ 50%, which request the ~ 10% changing of the basic 
parameters - average inter-instanton distance R and average instanton size p to restore 
chiral limit value of the quark condensate < qq > m =o and other important quantities as 
f n and m n to their phenomenological values. This work in the progress. 

In general, it was demonstrated, that instanton vacuum model is well working tool 
also beyond chiral limit and satisfy chiral perturbation theory. 
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